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Spatial 4a x 4a modulations, with a the lattice constant of Cu02 planes, or the so called 
checkerboards can arise from double-Q spin density wave (SDW) with Qi = (diTr/a, ±37r/4a) and 
Q2 = (it37r/4a, =l=7r/a). When multi-Q pair density wave, that is, the condensation of (i7-wave 
Cooper pairs with zero total momenta, it2Qi, ±2Q2, di4Qi, di4Q2, and so on is induced by the 
SDW, gaps can have fine structures similar to those of the so called zero-temperature pseudogaps. 
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High critical-temperature (high- Tc) superconductivity 
in cuprate oxides is an important and long standing issue 
since its discovery 1]. Many unconventional properties 
are observed in addition to high Tc'. the so called spin 
gap or pseudogap above H i, S H S H H 
4a X 4a spatial modulations called checkerboards, with a 
the lattice constant of Cu02 planes, of local density of 
states (LDOS) around vortex cores above [13, 
and below Tc |l5 0,^3 , the so called zero-temperature 
pseudogaps (ZTPG) [l^, and so on. The issue cannot be 
settled unless not only high Tc but also unconventional 
properties are explained within a theoretical framework. 

Since cuprates are highly anisotropic, thermal critical 
superconducting (SC) fluctuations, which would be di- 
vergent at nonzero Tc in two dimensions must play 
a role in the opening o f ps eudogaps. This issue is exam- 
ined in another paper |l8j . 

The period of checkerboard modulations is indepen- 
dent of energies. Their amplitude depends on energies; 
it is only large in the gap region. When the modulating 
part is divided into symmetric and asymmetric ones with 
respect to the chemical potential, the symmetric one is 
larger than the asymmetric one |l^. Fine structures are 
observed in ZTPG It is difficult to explain these ob- 
servations in terms of charge density wave (CDW). Sev- 
eral possible mechanisms have been proposed: Fermi sur- 
face nestin g |19l|. valence-bond solids E^, , pair den- 
sity waves I22I l23j , hole-pair [l^, 0, HBl oi" single- hole 
gner solids, and Wigner super solids 261 l27l. 
The purpose of this Letter is to propose another mech- 
anism: The spatial modulation of LDOS is due to spin 
density wave (SDW) and ZTPG is due to the coexistence 
of SDW and superconductivity or pair density waves in- 
duced by SDW. 

Cuprates with no dopings are Mott-Hubbard insula- 
tors, which exhibit antiferromagnetism. As dopings are 
increased, they exhibit the Mott-Hubbard transition or 
crossover and they becomes metals. High- Tc supercon- 
ductivity occurs in such a metallic phase. According 
to a previous theory which is consistent with a 

combined one of Hubbard's 29] and Gutzwiller's [s^ 
ones, a three-peak structure appears in the density of 



states; the so called Gutzwiller's quasiparticle band be- 
tween the lower and upper Hubbard bands (LHB and 
UHB). This is confirmed by the single-site approxima- 
tion (SSA) 31, 32] that includes all the single-site terms 
or by the dynamical mean- filed theory (DMFT) js^ . The 
three-peak structure corresponds to the so called Kondo 
peak between two subpeaks in the Anderson model or 
in the Kondo problem. A Kondo-lattice theory is use- 
ful to treat strong correlations in the vicinity of the 
Mott-Hubbard transition [s^. The superexchange in- 
teraction, which arises from the virtual exchange of pair 
excitations of electrons across LHB and UHB even in 
the metallic phase, favors an antiferromagnetic (AF) or- 
dering of (±7r/a, ±7r/a). Another exchange interaction 
arises from that of Gutzwiller's quasiparticles. Since 
Fermi surface nesting plays a role, this exchange interac- 
tion favors AF orderings of nesting wave numbers. Then, 
one of the most plausible scenarios for the modulating 
LDOS is that double-Q SDW with Q = (±7r/a, dz37r/4a) 
and (±37r/4a, ±7r/a) are stabilized and the 2Q modu- 
lation of LDOS is one of their second-harmonic effects; 
2Q = (±37r/2a, ±27r/a) or (±27r/a, ±37r/2a) is equiva- 
lent to (±7r/2a,0) or (0,±7r/2a). On the other hand, a 
magnetic exchange interaction can be a Cooper-pair in- 
teraction • According to previous papers 36, 37], the 
superexchange interaction as large as |J| = 0.1-0.15 eV is 
strong enough to reproduce observed SC Tc; theoretical 
Tc of (i7 wave are definitely much higher than those of 
other waves. Since the two intersite exchange interactions 
are responsible for not only antiferromagnetism but also 
superconductivity, the coexistence of antiferromagnetism 
and superconductivity or the competition between them 
must be a key for explaining unconventional properties 
of cuprate oxide superconductors. 

In order to demonstrate the essence of the mech- 
anism, we consider first a mean-field Hamiltonian 
on the square lattice, that is, non-interacting elec- 
trons in the presence of AF and (i7-wave SC fields: 
^ = Eka^(k)4a«ka + Haf + Hsc- The first 
term describes non-interacting electrons; T^(k) = 
—2t [(zos{kxa) (zo^{kya)\—At' (zo^ikxo) cos{kya) — fi^ with 
t and f transfer integrals between nearest and next- 
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nearest neighbors and fi the chemical potential, is the 
dispersion relation of electrons. We assume that t' /t = 
—0.3. The second term describes AF fields with wave 
number Q = (±37r/4a, ±7r/a) or (±7r/a, ±37r/4a): 

\<.aa' ^—x,y,z 

+ ^?4-Qa«ka') , (1) 

with a^^' the Pauli matrixes. A single-Q structure or 
the so called stripe is assumed for the sake of simplicity. 
The origin of real space is chosen in such a way that 
and A| are real and positive; the external filed is 
A^(R^) = 2A^ cos(Q • R^). The Brillouin zone is folded 
by periodic AF fields. When we take its origin at a zone 
boundary of a folded zo ne, electron pairs of k + /Q and 
— k + /Q can be bound j^^, with / being an integer. We 
assume the following dj-waye SC fields: 

k I 

+A^*a_k+^Qiak+^QT) , (2) 

with rjd^Ck) = cos{kxa) — cos{kya). The global phase 
of single-particle wave functions can be chosen in such a 
way that Aq and Aq are real and positive. We assume 
A^ = A_^ for / 7^ for the sake of simplicity, although 
we have no argument that other cases can be excluded. 
The homogeneous part of LDOS per spin is given by 

/5o(£) = - ^ E [^-(^ + ^^'k; 2^Q)],=o ' (3) 

kcr 

with 7/|t| +0, where Ga{£ + ^7,k; 2/Q) is the analyt- 
ical continuation from the upper half plane of 

G^iiSn^k; 21Q) = - Jdre''-^ (^rak-iQa(r)a^^Q^) , 

(4) 

with /3 = 1/kBT; we assume T = K so that /3 +oo. 
The modulating part with wave number 2/Q is given by 

kcr 

^g-i2iQ.R.^^(^ + i7, k; -2/Q)] , (5) 

Since A^ and Aq are real and Ai = A_/ for / 7^ 0, Gcr(^ + 
Z7, k; 2/Q) = Gaie^ij, k; -2/Q). Then, Eq. © becomes 
simple in such a way that 

P2iQ{e;'Ri) = 2cos(2/Q • Ri)p2/Q(e), (6) 

with 

P2iQ{e)=-^Y.^mG,{e + n,-k;±2lQ). (7) 

kcr 



The modulating part with (2/ + l)Q vanishes because the 
up and down spin components cancel each other. 

We assume that A^ = for |/| > 3 and we take 
a 5 X 2 X 2- wave approximation; we consider couplings 
among single-particle excitations such as electronic ones 
of «k,±. |0>, ^k±Q, ±cr ^k±2Q, ±cr |0) hole-like 

ones of a_k,±a |0), a_k±Q,±a |0), and a_k±2Q,±a |0), with 
|0) being the Fermi vacuum. Matrixes to be diagonal- 
ized are 20 x 20. The transformation diagonalizing them 
is nothing but a generalized Bogoliubov transformation. 
For the sake of numerical processes, nonzero 7 is as- 
sumed: We assume 7/|t| = 0.3 instead of j/\t\ +0. 

Figure Q show po(^) and P2q(^) in the presence of AF 
fields along the x axis and no SC fields; results do not 
depend on the direction of AF fields. A gap minimum 
is close to the chemical potential for = —1. This 

implies that the nesting wave number must be very close 
to Q for fi/\t\ = —1. We assume jii/\t\ = —1 in the 
following part. No fine structure can be seen in the low- 
energy part of po{e). The symmetric part of P2q(^) is 
much larger than the asymmetric one. The amplitude of 
CDW is small, because positive and negative parts below 
the chemical potential {e < 0) cancel largely each other. 

Figure [2I show po(^) in the presence of SC fields and 
no AF fields. In the absence of AF fields, even if A^ 7^ 
for / 7^ 0, there is no modulating part in LDOS. 

Figure [2I show po(^) and P2q(^) in the presence of AF 
fields along the x axis and SC fields; results do not depend 
on the direction of AF fields either. Gaps can have fine 
structures. The modulating part P2q(^) can be quite 
different between Figs. Q and [SI or in the the absence and 
presence of SC fields. 

In order to explain checkerboards and ZTPG in 
cuprates, various extensions have to be made. First of 
all, strong electron correlations in the vicinity of the 
Mott-Hubbard transition should be seriously considered; 
we had better use the so called d-p model or the t-J 
model . A Kondo-lattice theory can treat such strong 
electron correlations Non-interacting electrons in 

this Letter correspond to Gutzwiller's quasiparticles; ob- 
served specific heat coefficients as large as 14 mJ/K^mol 
imply \t\ 0.04 eV. External AF and SC fields in this 
Letter correspond to conventional Weiss's static mean 
fields due to the superexchange interaction and the ex- 
change interaction arising from the virtual exchange of 
pair excitations of Gutzwiller's quasiparticles. 

In the orthorhombic or square lattice, Qi = (±7r/a, 
±37r/4a) and Q2 = (±37r/4a, ±7r/a) are equivalent to 
each other. Since the Fermi surface nesting is sharp, 
double-Q SDW must be stabilized in cuprates rather 
than single-Q SDW; magnetizations of the two waves 
must be orthogonal to each other . We propose that 
checkerboards in the absence of SC order parameters are 
due to the second-harmonic effect of double-Q SDW. 

In the presence of double-Q SDW, superconductivity 
should be extended to include Cooper pairs whose total 
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FIG. 2: po(£) in the presence of SC fields (Ao / and 
= for / 7^ 0) and no AF fields; /i/|t| = —1 is assumed for 
the chemical potential. Since we assume nonzero ^ /\t\ = 0.3, 
po{e — 0) is nonzero. 



momenta are zero, ±2Qi, ±2Q2, ±4Qi, ±4Q2, and so 
on. Not only checkerboards but also ZTPG can arise 
in the coexistence phase of double-Q SDW and multi-Q 
superconductivity. The solid line in Fig. |3Ja)-(ii) resem- 
bles to observed fine structures of ZTPG, though single- 
Q SDW is assumed there. The observed ZTPG phase 
may be characterized by a precise comparison of po(^) 
and P2q(^) between observations and theoretical results 
for various parameters for double-Q SDW and multi-Q 
superconductivity. 

Although AF fluctuations are well developed in the 
checkerboard and the ZTPG phases, it is not certain 
that AF moments are actually present there. Checker- 
boards and ZTPG are observed by scanning tunnelling 
microscopy or spectroscopy, which can mainly see LDOS 
of the topmost Cu02 layer on a cleaved surface. A possi- 
ble explantation is that AF moments appear only in few 
surface Cu02 layers because of surface defects or disor- 



der. It is likely that disorder enhances magnetism. Ex- 
perimentally, for example, the doping of Zn ions enhances 
magnetism in the vicinity of Zn ions |4a; theoretically, 
magnetism is also enhanced by disorder '34^. 

It is reasonable that checkerboards appear around vor- 
tex cores because AF moments are induced there |lH : 
vortex cores can play a role of impurities for quasiparti- 
cles as doped Zn ions do. It is also reasonable that al- 
most homogeneous checkerboards appear in under doped 
cuprates where superconductivity disappears 0, [l5| : 
almost homogeneously AF moments presumably exist 
there. Inhomogeneous checkerboards can appear in un- 
der doped cuprates with rather low SC Tc's, when inho- 
mogeneous AF moments are induced by disorder. 

Even if the 2Q modulation in LDOS, P2q(£), is large, 
the 2Q electron density modulation or the amplitude 
of CDW caused by it is small because of the cancella- 
tion between positive and negative parts of P2q(^) below 
the chemical potential and the small spectral weight of 
Gutzwiller's quasiparticles, which is as small as with 5 
doping concentrations measured from the half filling, ac- 
cording to Gutzwiller's theory [s^. On the other hand, 
CDW can be induced by another mechanism in Kondo 
lattices 41j. In the vicinity of the Mott-Hubbard transi- 
tion, local quantum spin fluctuations mainly quench mag- 
netic moments; this quenching is nothing but the Kondo 
effect. The energy scale of local quantum spin fluctu- 
ations, which is called the Kondo temperature ksTx^ 
depends on the electron filling in such a way that it 
is smaller when the filling is closer to the half filling; 
ksTx — \St\ according to Gutzwiller's theory 30]. Then, 
CDW is induced in such a way that ksTx are smaller or 
local electron densities are closer to the half filling where 
AF moments are larger. Even if the amplitude of CDW 
is not vanishingly small in an observation, the observa- 
tion cannot contradict the mechanism of checkerboards 
and ZTPG proposed in this Letter. 

In conclusion, double-Q SDW with Qi = (±7r/a. 
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FIG. 3: (a) po{e) and (b) P2q(^) in the presence of AF fields of Ax/\t\ =0.3 and various SC fields; /i/|t| = —1 is assumed for 
the chemical potential. In (i), solid, dotted, and dashed lines are for (Ao, Ai, A2)/|t| = (1,0,0), (0,1,0), and (0,0,1). In (ii), 
they are for (Aq, Ai, A2)/|t| = (1, 0.5, 0), (1, —0.5, 0), and (1, 0, ±0.5), respectively; the result for (1, 0, —0.5) is the same as that 
for (1,0,0.5) within numerical accuracy. In (iii), solid and dotted lines are for (Ao, Ai, A2)/|t| = (l,dz0.5i,0) and (l,0,±0.5i), 
respectively; the result for (1, — 0.5i, 0) is the same as that for (1, 0.5z, 0) within numerical accuracy, and so on. 



±37r/4a) and Q2 = (±37r/4a, ±7r/a) must be responsible 
for the so called 4a x 4a checkerboards in cuprate ox- 
ide superconductors, with a the lattice constant of Cu02 
planes. Not only Cooper pairs with zero total momenta 
but also those with ±2Qi, ±2Q2, ±4Qi, ±4Q2, and 
so on are possible in the SDW phase. The so called zero 
temperature pseudogap phase must be a coexisting phase 
of the double-Q SDW and the multi-Q condensation of 
dj-waye Cooper pairs. 

The author is thankful for discussion to M. Ido, M. 
Oda, and N. Momono. 



* Electronic address: fohkawa@phys.sci.hokudai.ac.jp 
[1] J.G. Bednortz and K.A. Miiller, Z. Phys. B 64, 189 
(1986). 

[2] H. Yasuoka et al.. Strong Correlation and Superconduc- 
tivity, Springer Series in Solid State Science Vol. 89 
(Springer- Verlag, Berlin, New York 1989), p. 254. 

[3] H. Ding et al.. Nature 382, 51 (1996). 

[4] J. M. Harris et al., Phys. Rev. B 54, R15665 (1996). 

[5] A. G. Loeser, et al.. Science 273, 325 (1996). 

[6] A. Ino et al., Phys. Rev. B 65, 094504 (2002). 

[7] Ch. Renner et al., Phys. Rev. Lett. 80, 149 (1998). 

[8] M. Oda et al., Physica C 281, 135 (1997). 

[9] T. Nakano et al., J. Phys. Soc. Jpn. 67, 2622 (1998). 
[10] T. Ekino et al., Phys. Rev. B 60, 6916 (1999). 
[11] J.E. Hoffman et al.. Science 295, 466 (2002). 
[12] M. Vershinin et al.. Science, 303, 1995 (2004). 
[13] T. Hanaguri et al.. Nature 430, 1001 (2004). 
[14] C. Howald et al., Phys. Rev. B 67, 014533 (2003). 
[15] N. Momono et al., J. Phys. Soc. Jpn. 74, 2400 (2005). 
[16] K. McElroy et al., Phys. Rev. Lett., 94, 197005 (2005). 
[17] N.D. Mermin and H. Wagner, Phys. Rev. Lett. 17, 1133 
(1966). 



[18] F.J. Ohkawa, |cond-mat/0508"344l 

[19] H.C. Fu et al., cond-mat/040300r 

[20] S. Sachdev, Science 288, 475 (2000). 

[21] M. Vojta, Phys. Rev. B 66, 104505 (2002). 

[22] D. Podolsky et al., Phys. Rev. B67, 094514 (2003). 

[23] H.-D. Chen et al., Phys. Rev. Lett. 93, 187002 (2004). 

[24] H.-D. Chen et al., cond-mat/040232. 

[25] H.-D. Chen e tjiL, Phys. Rev. Lett. 89, 137004 (2002). 

[26] Z. Tesanovic, cond-mat/0405235 

[27] P. W. Anderson, cond-mat/0406038l 

[28] F.J. Ohkawa, J. Phys. Soc. Jpn. 58, 4156 (1989). 

[29] J. Hubbard, Proc. Roy. Soc. London Ser. A 276, 238 

(1963); A 281, 401 (1964). 
[30] M.C. Gutzwiller, Phys. Rev. Lett. 10, 159 (1963); Phys. 

Rev. A 134, 293 (1963); A 137, 1726 (1965). 
[31] F.J. Ohkawa, Phys. Rev. B 44, 6812 (1991). 
[32] F.J. Ohkawa, J. Phys. Soc. Jpn. 60, 3218 (1991); 61, 

1615 (1992). 

[33] A. Georges an d G. Kothar, Phys. R ev. B 45, 6479 (1992). 

[34] F.J. Ohkawa, |co nd-mat/0510377| To be published in J. 
Phys. Soc. Jpn. 

[35] J.E. Hirsch, Phys. Rev. Lett. 54, 1317 (1985). 

[36] F.J. Ohkawa, Jpn. J. Appl. Phys. 26, L652 (1987); J. 
Phys. Soc. Jpn. 56, 2267 (1987). 

[37] F.J. Ohkawa, Phys. Rev. B. 59, 8930 (1990). 

[38] A superconducting state with / for / / is noth- 
ing but a pair- density wave state or a generalized Fulde- 
Ferrel-Larkin-Ovchinnikov (FFLO) state induced by AF 
fields. Electron pairs of k+ (/ + |)Q and -k+ (/ + |)Q, 
whose total momenta are (2/ + 1)Q, are also possible. 
However, they are ignored here; they should be consid- 
ered in case of CDW. 

[39] E.G. Zhang and T.M. Rice, Phys. Rev. B 37, R3759 
(1988). 

[40] J.W. Loram, et al., Phys. Rev. Lett. 71, 1740 (1993). 
[41] F.J. Ohkawa, J. Phys. Soc. Jpn. 67, 535 (1998); Phys. 

Rev. B 66, 014408 (2002). 
[42] H. Alloul et al., Phys. Rev. Lett. 67 (1991) 3140. 



